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O : 

The effect of random surface roughness on hydrodynamics of viscous incompressible liquid is 
. discussed. Roughness-driven contributions to hydrodynamic flows, energy dissipation, and friction 

force are calculated in a wide range of parameters. When the hydrodynamic decay length (the viscous 
wave penetration depth) is larger than the size of random surface inhomogeneities, it is possible to 
replace a random rough surface by effective stick-slip boundary conditions on a flat surface with two 
constants: the stick-slip length and the renormalization of viscosity near the boundary. The stick- 
slip length and the renormalization coefficient are expressed explicitly via the correlation function 
of random surface inhomogeneities. The effective stick-slip length is always negative signifying 
the effective slow-down of the hydrodynamic flows by the rough surface (stick rather than slip 
motion). A simple hydrodynamic model is presented as an illustration of these general hydrodynamic 
results. The effective boundary parameters are analyzed numerically for Gaussian, power-law and 
exponentially decaying correlators with various indices. The maximum on the frequency dependence 
ryj ' of the dissipation allows one to extract the correlation radius (characteristic size) of the surface 

inhomogeneities directly from, for example, experiments with torsional quartz oscillators. 

PACS numbers: 47.10.+g, 68.08.-p, 46.65.+g, 81.40.Pq 
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O I I. INTRODUCTION 

Progress in micro- and nanotechnology requires better understanding of boundary effects. For hydrodynamic 
^ ■ microflows, this means better understanding of stick-slip motion near solid walls and, in particular, information on 
the dependence of the slip (or stick) length on the properties of the walls. Despite the fact that similar issues were first 
raised more than a hundred years ago [Q, ^, 0], the slip length remains one of the least known transport coefficients. 
Traditionally, the most detailed information on the boundary slip is available for rarefied classical gases [ 
applications to vacuum technology, high altitude flights, and space research. More recently ||, ^, [h], |lT|, |i"2, 
liquid 3 He has become an important source of information on surface slip. This is not surprising since, in contrast 
to classical gases, one can easily vary the quasiparticle mean free path in 3 He by changing temperature thus allowing 
■ experiments in a wide range of Knudsen numbers. 

The conventional theory of boundary slip assumes that the slip length is proportional to the bulk mean free path 
H £{,, C s i — aCb, and ignores small-scale surface inhomogeneities. Obviously, this approximation is too crude. The 
hydrodynamic flows near the walls strongly depend on geometry of surface inhomogeneities |l6|]. Recent analysis of 
slip near a model surface with periodic irregularities demonstrated [l7| that the effective slip length C e f f contains not 
only the bulk component aCb but also the contribution from the averaged surface curvature R, £~ff — or 1 C^ 1 — Rr 1 , 
An application of the corresponding boundary condition to several types of curved walls |Q resulted in an interesting 
expression for an effective slip length which could, under certain circumstances, be equivalent to large-scale surface 
roughness. However, the results jl7], [l8| were obtained for few special types of regular surface inhomogeneities only. 
In the case of micro- and nanoscale defects, it is more realistically to suggest that surfaces have random corrugation. 
What is more, in some cases, especially in the hydrodynamic limit Cb — ► 0, it is not clear how to use the effective 
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boundary parameters of Refs. 17, [L8|. 

Below we derive an effective stick-slip boundary condition which would reproduce hydrodynamic flows with Cb = 
near rough walls with small-scale random inhomogeneities. Since the hydrodynamic calculations near inhomogeneous 
walls are extremely complicated | fl6| , it is highly desirable to map this problem onto the system with simple flat 
surface geometry with some effective boundary condition. This boundary condition should contain information about 
geometrical and statistical properties of the real corrugated surface and ensure a proper behavior of hydrodynamic 
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FIG. 1: General geometry of the model. 



variables. The derivation of this simple boundary condition is the main goal of the paper. Below we show that this 
boundary condition contains two effective parameters: the effective stick-slip length and renormalized viscosity. We 
also demonstrate that the results for attenuation in torsional oscillator experiments can provide valuable information 
about the statistical type of surface inhomogeneities and give the values of the main geometrical parameters of surface 
roughness. 

In the next Section, we present the main hydrodynamic equations and find the stream function in systems with 
random rough walls (details of the derivation are given in Appendix A). Comparison of the hydrodynamic results 
with those for the stick-slip motion allow us to get the expression for the effective stick-slip parameters in Section III. 
For clarification of the physical meaning of the parameters in the somewhat unexpected effective boundary condition, 
we present a simple hydrodynamic model for a boundary layer in Appendix B. Section IV contains analytical and 
numerical results for surfaces with various statistical types of inhomogeneities. Summary and conclusions are presented 
in Section V. 



II. HYDRODYNAMIC FLOWS ALONG ROUGH WALLS 



To determine an effective slip/stick length, one has to solve an appropriate hydrodynamic problem with a boundary 
condition on a random rough wall and to compare the results with those for a similar problem with a slip boundary 
condition on a smooth wall. Several "typical" hydrodynamic problems ||l9f have been generalized recently in order 
to cover boundaries with slight roug hness @, [HJ || §| ||. For our purposes, the most appropriate problem is the 
problem of hydrodynamic flows excited by tangential oscillations of a rough wall. The advantages are the convenience 
of the experimental setup with a standard transverse oscillator, a choice of several observables such as hydrodynamic 
velocity and two components of the shear impedance, and the presence of an extra variable - frequency uj - that 
allows one to vary the ratio of the hydrodynamic decay length to the size of wall inhomogeneities. Since this problem 
has already been studied in Ref. p0[ , though by a different method, we will only briefly outline our hydrodynamic 
formalism in Appendix A and present some additional results. 

We consider semi-infinite viscous fluid restricted by a rough solid wall. For simplicity, roughness is assumed one- 
dimensional with profile described by a random function Y = E (X) with the zero mean value, (H (X)) = 0. The wall 
is homogeneous in Z-direction (see Fig. This inhomogeneous surface is characterized by two length parameters: 
the average amplitude h and correlation radius (size) R of surface inhomogeneities. We consider the case of slight 
roughness, 

e = h/R<^l. (1) 

In other situations, any general description of hydrodynamic flows near rough walls is virtually impossible. 
The wall oscillates in ^-direction with the velocity 



U (t) = e^uo cos (ujt) (2) 

The hydrodynamic boundary condition is the condition of zero velocity V on the wall in the reference frame in which 
the wall is at rest: 



V^X- J U (t) dt,Y = 3 (x - J U(t)dt^j = 0, 



(3) 
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Two important hydrodynamic length scales are the decay length (or the viscous wave penetration depth), 8, and 
the amplitude of the boundary oscillations, a, 

5 = \j2i//u), a — uq/uj, (4) 

where v = rj/ p is the kinematic viscosity. 

It is convenient to choose h, R and the amplitude of the wall velocity uq as the scaling parameters and introduce 
dimensionlcss variables as 



v = V/uo, x = X/R, y = Y/R, £(x) = E(X/R)/h. (5) 

When the fluctuations of are statistically independent and the higher momenta can be expressed through the 
second one, the random surface roughness is actually described not by the unknown random function £(x) with the 
zero average, but by the correlation function £ (x): 

C(x) = {£(x 1 )S(x 1 + x)) = ^J Z(xi)t( Xl + x)dxi, (6) 

«(*xK(0> = 27rf(k + *4K(**), 

where A is dimensionless flat surface area of the wall. Experimentally the correlation functions £ (x) [or its Fourier 
image, also called the power spectrum, £ (k x )] can exhibit different types of long-range behavior and assume various 
forms p5fl . Particular examples of the surface correlators are analyzed in Section IV. Note, that in our dimensionless 
notations (^|) the correlation radius of the surface inhomogeneities is equal to 1. 

The liquid is considered incompressible, divv =0. In variables (0), the dimensionless Navier-Stokes equation can 
be written as 

1 d rot v 9 m 

-V 2 rotv = 3?[(rotvV)v-(vV)rotv] (7) 



cjo dt 

where the characteristic frequency u>o and the Reynolds number -ft are 

-o = ^=^|- (8) 
R l v Rujq 

(the inverse frequency parameter uj^ 1 is often called the diffusion time of vorticity) Since the first term in Eq. (Q) 
has an order of (w/wq) K>t v, the hydrodynamic flows are characterized by the dimensionless parameter 



A = ^/uj/ujo = V2R/S (9) 

which describes the ratio of the size of inhomogeneities R to the hydrodynamic decay length i5. Two dimensionless 
parameters, e and A, are the main parameters of the problem. 

Below we consider the linearized Navier-Stokes equation without the nonlinear term in r.h.s. of Eq. (Q). For small 
frequencies, lo/luq <C 1, this linearization is justified for small Reynolds numbers 5ft -C 1. In the opposite limit of high 
frequencies, lu/ujq 3> 1, this requires smallness of the amplitude of oscillations a in comparison with the tangential 
size of surface inhomogeneities R at arbitrary Reynolds numbers 5R, a/R <§; 1 |l9| ]. The linearized Eq. (0) for rotv 
can be, as usual, rewritten as the fourth order differential equation for the scalar stream function ip (x,y), 



dip dip 
dy ' v dx 



- = (io) 



In our problem, all hydrodynamic variables contain harmonic time dependence. After the transformation to the 
coordinate frame oscillating with the wall, the hydrodynamic equations and boundary conditions for the stream 
function acquire the form 

-zA 2 VV - VV = 0, (11) 
dip{x,et,{x)) _ di/j(x,e£(x)) _ 

8y " ' dx _U ' [ ' 

ip(x,oc) — const. (13) 
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The solution of the linearized Navier-Stokes equation JIH ) — fll5| ) is quite difficult because the boundary condition 
( |l2|) involves the rough wall with random inhomogeneities. Using a coordinate transformation y — > y — £ (x) , we can 
reduce the Navier-Stokes equation to an equivalent equation with the boundary condition on the perfect flat wall. 
However, this new equation, as a result of the transformation-driven change in derivatives, acquires several additional 
terms Vip that involve the combinations of derivatives of ip and the random function £ (x). To deal with these terms, 
we find the explicit form of the Green's function with the proper boundary condition. Then the problem reduces to 
a rather transparent integral equation 



''(/,,.,/) = ,-,:„,,(/,■.,/)+ / <h,'G (/.-.,..//.//) / ^V(k x -k' x ,y')i>(k' x ,y') 



(14) 



This procedure and the explicit expressions for the unperturbed inhomogeneous solution tpinh {k x , y) , the perturbation 
V, and the Green's function are given in Appendix A. In some sense, we shifted the difficulty from the boundary 
condition to the bulk equations with random sources of the special form. Note, that Eq. (|lj) is still exact and, in 
principle, could be solved without the perturbation theory. The explicit form of Green's function is such that one 
can extract the main part of the solution in the closed form. Another possible approach to Eq. ( |l4| ) is to apply the 
Wiener- Hermite functional expansion pfl, E7|. 



Here we solve Eq. ( |14| ) by iterations as an expansion in the small parameter: ip = ipo + eV^i + £ ip2 
three terms for the stream function have the following forms 

2tt 



■00 (k x ,y) = -r-<5 (fc x )exp(iAy) , 

iX 



i\y 



iX 



S 2 ~ Si 



(^(k x ,y)) = S(k x ) / dk' x ((k' x ) 



Sl + s 2 
iX 



(iXe iXy + Sl e SlV - s 2 e S2V ) 



The first 

(15) 
(16) 

(17) 



where we exclude uninteresting constant terms and 



si = - \k x \ , s 2 = v&l - = -Oi + ifi, 
X = e"/ 4 A. 



(18) 



Since for further calculations we need only the expression for i\j 2 which is averaged over the random surface inhomo- 
geneities, Eq. ( |l7| ) gives only the compact expression for (k x , y))- 

These expressions for the stream function provide the roughness-driven corrections for the velocity and rate of 
energy dissipation (see Appendix A): 



(v x ) = Re | e «(A»-«0 [j + iAe 2 £i ] j ; {Vy) = Qj 



h = 



^C(k x ){ Sl+S2 -tX/2} 



g-^[l + e>A^] 



£ 2 = [ -C(tA)<f>(t), $(t) = l-\fjt+t±-t*. 



(19) 
(20) 

(21) 
(22) 



The equation for the energy dissipation is averaged over both the surface roughness and the period of oscillations. 
This expression is similar to the result of Ref. [g0| obtained with the help of the Rayleigh perturbation method. 

Stream function also allows one to find corrections to the roughness-driven friction force. These calculations should 
be done more carefully than for standard flat geometry: the friction force is parallel to the actual surface and, in the 
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case of the oscillating rough wall, has both components F x and F y . One should also take into account the y-component 
of velocity, which is absent in the case of flat geometry. Straightforward calculation for the averaged square of absolute 
value of dimensionless friction force give 



' R 



A 2 



e 2 A 2 



dtC,{kt)(j> 2 it) 



(23) 



This expression is different from a simple experimental definition of the effective friction force F e tt = —Q/uq 
At low frequencies (large decay lengths, A< 1), Eqs. (|l9|) , (^lj) expressions for parameters reduce to 



OO 

/dk x _ 

o 



0(A), 



£ 2 = 1 + (A hi A), 
and the equations for the velocity and attenuation acquire the following form: 



(v x ) = Re 



a i(\y-u>t) 



(oo 
2 J ^((k x )k x + 



(A) 



Q=-ffA[i + A V(l + 0(AlnA))] 



(24) 
(25) 

(26) 
(27) 



The fact that the main term in £2 is equal to 1 is due to our choice of the normalization of the correlation function 
in Eq. (|j) as £ (x = 0) = 1 (see also Section IV). 
In the opposite limit of high frequencies A>1, 



t]Uq A 
2i? V2 



dk 



Qo 



(28) 



This result has a simple physical explanation. In this limit, the decay length is much smaller than correlation radius 
(size) of the wall inhomogeneities R. As a result, the dissipation occurs in a very narrow layer near the wall within 
which the wall can be considered as almost flat. Then the correction to dissipation stems simply from the increase in 
the surface area relatively to the flat boundary 



Q „ -2^41 l dA = JL I / y/T+^dX 



2R y/2L 2 



2R y/2L 



(29) 



Equation fl28J) is simply the combination of the first two terms in the Taylor expansion of Eq. (|29|) in small e. 

In principle, it is possible to slightly modify our problem by considering a torsional quartz crystal oscillator with 
density p s , thickness d. If such a resonator has a rough solid- fluid interface, the frequency shift Slo of the resonance 
frequency O acquires an additional roughness-driven component which can be described within the above formalism 
and should be given by the similar equations. Such a frequency shift for a transverse oscillator is [^(J 



V2Rp s d 



—C(k x ) 



^k x + A 4 + kl - A + V2k x 



We do not want to dwell on this issue; our interest in focused mainly on the roughness-driven corrections to the 
hydrodynamic flows and dissipation. 



III. EFFECTIVE STICK-SLIP BOUNDARY CONDITIONS 



The main aim of the paper is to find when and to what extent flows near random rough surface are equivalent to 
stick-slip motion with some effective stick-slip boundary conditions near flat surfaces, 



Re { v x (x, 0, t) - % ^L%M \ = Re («"*-) 



R 



Oy 



(30) 
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where the effective stick-siip iength C e ff, in order to simplify the applications of the results, is introduced with the 
proper dimensionality of length while all other variables are still dimensionless, Eq. (|J). With this boundary condition 
on a flat wall, the velocity field is 



Vx {y, t) = Re 



(31) 



Since the roughness-generated corrections for velocity are small, the comparison between Eq.(|3l|) and Eqs.([T^) , ( |26| ) 
is possible only when A£ e ff/R <C 1, i.e., only for relatively large decay lengths (low frequencies), 



v x {y, t) ~ Re 



1 + e^^KCeff/R 



(32) 



In this case, the comparison with the roughness-driven correction for the velocity at low frequencies, Eq. (|26|), yields 
the following simple expression for the effective stick-slip length C e ff = Re 2 li: 



R 



oo 

/ C ij^x) k x 

J 7T 



(33) 



means that the rough boundary roughness causes effective slow-down of the liquid, i.e., 
is the stick length rather than the slip length. In other words, there is an additional 



The negative sign in Eq. ([ 
that the coefficient £ e ff (| 
roughness-induced friction. 

The condition AC e ff/R 1, when one can replace the rough wall by a stick-slip boundary condition is equivalent 

to 



AC 



eff 



R 



Ae 2 



h2 , 

R6 <<:1 



(34) 



Surprisingly, the effective boundary condition ( |30| ) , (|3^) , taken by itself, cannot emulate the roughness-driven atten- 
uation (pl|). The reason is the presence normal flows near the boundary, v y (x), which are completely absent within 
the effective stick-slip description ( |30| ) , ( |3l| ) in which v y — 0. The attempts to modify the boundary condition ( |30| ) 
so that to reproduce both the velocity and attenuation correctly by, for example, introducing a two-component or 
complex stick-slip length, fail. In order to emulate the correct behavior of liquid near a rough wall, one has not only 
to introduce the stick-slip length (E^Ol) , (|33|) , but also to renormalize the viscosity near the wall as 



Vefj (y) = v [i + PS (»)] , 



where renormalization parameter (3 is equal at small A to 







R 



71^ 



(35) 



(36) 



The effective boundary conditions fl30|) , (|33|) , (|35|),(p6[) are the main result of this paper. These conditions allow one 
to replace the random rough boundary by an equivalent problem with the flat boundary and the effective stick-slip 
length and renormalized viscosity. The necessity of the renormalization of the viscosity means that the rough surface 
slows the flow down and changes the attenuation. Usually, the slip boundary condition is understood in terms of the 
existence of a peculiar thin slip boundary layer with the thickness of the order of the mean free path and with the 
properties that are somewhat different from the rest of the liquid. In the case of the rough walls, one should not 
only introduce the effective stick-slip layer with the thickness that is determined by S and R, but also to renormalize 
the viscosity in this layer explicitly. A simple physical model that clarifies the meaning of the effective parameters is 
given in Appendix B. 



IV. COMPARISON FOR DIFFERENT TYPES OF RANDOM INHOMOGENEITIES 



In this Section we address the question whether it is possible to extract information on the properties of the rough 
surface from the frequency dependence of attenuation of transverse oscillations. Statistical properties of the random 
surface are described by the correlation function of surface inhomogeneities, S (X) — h 2 ((x), x = X/R, Eq. (^J) . 
Experimentally, the correlation function can exhibit different types of long-range behavior and can assume various 
forms fl25l. 
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TABLE I: The position of the maximum of the funcion ( |37| ) and the value of stick-slip length Eq. ( |43[ ) for different types of 
the surface correlation function. 



J 

8 


Correlator type 


Form, £ (x) 


1. 


Gaussian 


1 2\ 

exp [ — x j 


2 


Power-law 


(l+x2) -(M+l/2) 


2a 


/i = 1/2: Lorentzian 


(l+x 2 )- 1 


2b 


fi — 3/2: Staras 


(l+x 2 )- 2 


3 


Power-law Fourier image 




3a 


v — 1/2: exponential 


exp (— \x ) 


3b 


v = 3/2 


(1 + x ) exp (- a; ) 



Fourier image, £ (fc) 



An 



-tx (A < 1) 



y^exp 

2f-ir( M +i/2) 



-fc 2 /4) 



7rexp (— |fc|) 
f (l + [fc[)exp(- 
20F^^ ( l +fc2 ) 

2(l + fc 2 )" 1 
4(l + fc 2 )- 2 



1*1) 

{-+1/2} 



1.293 

1.320 
1.825 

no max 
1.238 



4/^F 
4 r( M +i) 
v^F r( M +i/2) 

2 
3 

2 r(./+i/2) i 
v^F r(i^) ^-1/2 

oo 

4/tt 



Three broad classes of the correlation functions C (^) an d their Fourier images C, (k) (the so-called power density 
spectral function, or power spectra) are summarized in Table [j]. For better comparison, all the correlators are 
normalized in the same way, £ (x = 0) = 1. Note, that this normalization differs from the one used in Ref. [^8) for 
conductivity of ultrathin films: the natural reference point for the conductivity was its value in the limit kR — > and 
all the correlation functions in Ref. pc| have been normalized using £ (k = 0) = 1. For the hydrodynamic problem in 
this paper, the normalization £ (a; = 0j = 1 provides a better reference. 

The most commonly used correlation function, namely, the Gaussian correlator, is listed first. The next class of the 
correlation functions covers power-law correlators with the exponentially decaying Fourier images (power spectra), 
|/c| M (|/c|). Here, the most widely used are the Lorentzian correlator (index /i = 1/2) and the Staras correlator 
(/i = 3/2). The third class of the correlation functions includes the conjugate correlators, namely, the exponentially 
decaying correlators with the power-law spectral function C(^)- I n our dimensionless notations, Eq. (^|), all the 
correlators have the correlation radius equal to one. 

The most convenient observable is the frequency dependence of the relative attenuation, Eq. ( pl| ) : 



oo 

Ae z Q J 7T 



(37) 



(f>(t) = l- y vT+t^-t 2 . 

In the limits i< 1 and t ^> 1, the function 4> (t) has the following asymptotic expansions: 



(38) 



•(*) 



i 2 /2-i 4 /8, 
1 - l/y/2t, 



t < 1 
t > 1. 



(39) 



Note, that the piecewise continues function, defined by the expressions in Eq. (^9|) connected at the point t — \/2, 
gives a good approximation for <j) (t) in the whole range of t. This can be useful in simple approximations of the 
integral (|37|). The integral ( |37| ) can be conveniently split into two parts, AFi and Ar2, which correspond to the 
contributions from small and large t. 

In the laminar limit, A <^ 1, the main contribution comes from large t : 



Ar - Ar 2 = A - A 2 In (A) fy/2 + O (A 2 ) 



(40) 



The first two terms in this expression are the same for the correlators of all types. Therefore, in the low-frequency 
limit with large decay length it is impossible to distinguish statistical properties of different surfaces. The physical 
reason is obvious: large-scale attenuation processes on the scale of decay length S are not very sensitive to the details 
of surface inhomogencitics with the size R S. 

The situation is different in the opposite case of large A. In this limit for Gaussian and power-law correlators with 
the exponential power spectra (types 1 and 2 in the Table |), the contribution from large t to the integral j3^ ) is 
exponentially small. An estimate of the contribution from small k yields 



Ar - Ari 



dk k 2 ((k) 
~ 2A 



1 

2A 



dx 2 



x=0 



a 

A' 



(41) 



where a = 1 for the Gaussian correlator and a = fi + 1/2 for the power-law correlators. 
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FIG. 2: Correction to the energy dissipation rate, Ar as a function of the frequency parameter A for Gaussian and ^-correlators 
in log-log scale. 



For the correlators with the power-law power spectrum (correlators of the type 3 in the Table |), the contribution 
from large t, AT 2, is 

Ar 2 ~ a [ 1/A t2v ~ ldt „ a 1 - 2 ^. 

Jo + 

The contribution from small t, ATi, strongly depends on the value of the exponent v. If < v < 1, then the value of 
Ari is determined by the upper limit of the corresponding part of the integral and it is also proportional to A 1 " 2 ". 
If v > 1, then the first terms in the Taylor expansion for tfi(t) yields a convergent integral proportional to A -1 , while 
the rest gives the terms with the smaller exponent K 1 ~ 2lJ : 



AT 1 



dk 



(l + fc 2 ) { " +1/2} 



fc 2 k 4 
2A ~ 8A 3 



1 of 1 



A VA 2 "- 1 



Thus, the energy dissipation rate for the correlators with the power-law power spectrum is determined by the value 
of the index v : 

Ar(A,„)^j A A 1 t-»<"< 1 - (42) 
A \ v > 1. 



Comparison of the asymptotic behavior of the function Ar (A) for small and large A, Eqs. ( j40|) — (|42|), indicates 
that this function should have a maximum at A — y/2R/S ~ 1 except for the correlators with small v. In experiment, 
the position of this maximum on the frequency dependence of the attenuation can become a direct measurement of 
the correlation radius (size) of the surface inhomogeneities R. 

The numerical results are summarized in Fig. which presents the functions Ar (A) for various correlators. Nu- 
merical values of the position of the maximum for Ar (A) for various correlation functions are presented in the 
Table I 

The last column in the Table describes the dimensionless roughness-driven stick-slip length l\ — £ e ff/e 2 R, Eq. 
(B3T), for various correlators at small A, 



2 r^c(k x )k x 

Jo * 



MA«1) = 2/ —C(k x )k x . (43) 
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V. 



SUMMARY 



In summary, we calculated roughness-driven contributions to the hydrodynamic flows, energy dissipation, and the 
friction force in a wide range of parameters. We also investigated the possibility of replacing a random rough surface 
by a set of effective stick-slip boundary conditions on a flat surface. Such a replacement is highly desirable for analysis 
of experimental data and/or simplification of hydrodynamic computations in microchannels. The replacement turned 
out to be possible when the hydrodynamic decay length (the viscous wave penetration depth) is larger than the size 
of random surface inhomogeneities, Eq. (|34|) . The effective boundary conditions contain two constants: the stick-slip 
length and the renormalization of viscosity near the boundary. The stick-slip length and the renormalization coefficient 
are expressed explicitly via the correlation function of surface inhomogeneities. The corresponding expressions are 
quite simple and can be easily used for analysis of experimental data or in hydrodynamic computations. The effective 
stick-slip length is negative meaning the effective average slow-down of the hydrodynamic flow by the rough surface 
(stick rather than slip motion). 

For better understanding of the results, in Appendix B below we present a simple hydrodynamic model that 
illustrates our general hydrodynamic calculations. 

In the process of the derivation of the effective boundary condition, we reduced the Navier-Stokes equation with the 
boundary condition on a random rough wall to the exactly equivalent closed integral equation with the homogeneous 
boundary condition on the ideal flat wall. All the information on the surface roughness is contained in the kernel of 
this integral equation. The equation can be solved by standard methods. 

The effective boundary parameters are analyzed numerically for three classes of surface correlators including the 
Gaussian, power-law and exponentially decaying correlators. The energy dissipation near the rough surface is calcu- 
lated as a function of frequency for these types of the correlation functions. The position of the maximum on the 
frequency dependence of the dissipation allows one to extract the correlation radius (characteristic size) of the surface 
inhomogeneities directly from, for example, experiments with torsional quartz oscillators. The stick-slip length is also 
evaluated numerically for all three classes of surface correlators. 

The next step should be the evaluation of the effective stick-slip length for ultrathin flow channels of the thickness L 
for which L is expected to gradually replace the decay length 8 in the expressions for the slip length. Another desirable 
development would be the incorporation into the effective slip length of both surface and bulk scattering processes 
beyond the simple Matthiessen's rule in the same spirit as recent calculations for helium flows in microchannels p9[. 



The work is supported by NSF grant DMR-0077266. 

APPENDIX A: SOLUTION OF THE NAVIER-STOKES EQUATION FOR FLUIDS RESTRICTED BY 

RANDOM ROUGH WALLS 

First, we assume that all the variables have the harmonic time dependence, exp (— iujt), transform the linearized 
Navier-Stokes equation (M) to the non-inertial coordinate frame in which the wall is at rest, 
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exp (— iut) dt, 



and introduce the stream function ip{x,y) as 



dip 



dip 
dx 



(Al) 





the boundary condition 



(I) 



can be rewritten as the 



(A2) 



dip (x,e£ (x)) 



= 0, 



(A3) 



dx 



(A4) 



10 



The difficulty in solving Eqs. (A2) — (A4) originates from the presence of a random function £ (x) in the boundary 
condition. The next step is the coordinate transformation 



x -> x, y -> y - e£ (x) , 
that flattens the wall, making the boundary condition ( [L2| ) simple, 



9y 



1. 



<9x 



e^x (x) , ^>(x, oo) = const. 



The change in derivatives introduces the additional term V (£,d x )ip into the r.h.s of Eq.(A2), 

-ik 2 V 2 i) - v 4 v = V (f , <9 X ) ^ 

where 



eVi + e 2 F 2 + e 3 l/ 3 + e 4 F 4 



r 

Vl^ = - (Z^Llpyx + y?ixxi>y + ^xxxlpyx + 6 txx^Pyxx + £, xxxx^Py ^^xxlpyyy ~t~ ^^x^Pyxxx ^^x^Pyyyx) 

v 2l p 



^ie ^yy "I - ^£,xxx£,x^Pyy H~~ ^-^^x^xx^Pyyx H~~ 6^ Ipyyxx ~h 2^ ^yyyy H - ^^a;a; V^y 
2£a; {^^xx^yyy H~~ ^^x^yyyx) j 
£e ^yyyyi 



(A5) 
(A6) 

(A7) 
(A8) 



and lower indices denote the differentiation of the functions £ and ^. 

The simplicity of boundary condit ion s i n ne w coordinates allows us to find the Green's function G(x — x',y, y') for 
the operator in the l.h.s. of Eqs. (|A2|) , (A7). With the help of this Green's function, our initial problem with a 
boundary condition on the random rough surface reduces to the compact integral equation: 



dy> G{k x ,y,y>) J (k x - k' x , y') $ (k' x , y') 

-oo 



(A9) 



where we performed the Fourier transformation in x-direction (in the new coordinate frame, the geometry of the 
boundary is independent of x) , and 



tyinh (k x ,y) 



dxe 



(A10) 



is a solution of Eq. (A2) with V — and boundary conditions (A3),(A4). With this definition of tpmh y) , 
the Green's function satisfies the homogeneous boundary conditions on the wall. Note, that Eq. ( A9) is an exact 
equivalent of our initial problem with the random rough wall and, in principle, can be solved for an arbitrary function 
£(x). 



The function ipi n h (k x , y) is determined by the characteristic equation for the operator in l.h.s. of Eq. ( A2) 

.s 4 - (2fc2 _ Z A 2 ) s 2 - (ik 2 x A 2 - fc 4 ) = 0. (All) 



This equation has four solutions: 



Si 



±si,±s 2 ; 



(A12) 



a,f3 = 



- | ^ | , s 2 = yfef - «A 2 = —a + i(3, 
-i= V / V^TA^±fc2>0. 



We are interested only in the functions t/jinh (k x ,y) that decrease at y — > oo. Therefore, the general solution of the 
homogeneous Eq.(A2) with the boundary condition (A3) has the form 



iPinh (kx,y) = ^5 (k x ) [e iX y 1] + 

l\ S2 — S\ 



S2C 



Sl e S2V ] 



(A13) 
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and contains the contribution without e, ipQ (k x ,y), and the term linear in e. Similar calculations yield the Green's 
function: 



G{k x ,y,y') 



1 



2ik 2 



1 

si 



i(v+v') 



(A14) 



iA 2 (s 2 



si 



si(y+y') _|_ gS2(a+2/') _ fi siJ/+s 2 i/ _ p siy'+s 2 y 



The last result can be also obtained by noticing that our Green's function is proportional to the difference between 
the Green's functions for the two-dimensional Laplace and Hclmholtz equations with the same boundary conditions: 



G(r,r') = \- 2 {G L -G H ), 



G L (r,r') = -^ln(R s /Ri), G 2 {v,v') = - (XR S ) - H^> (XRj 

Z7T 4 L 



Rsj = \]{x - x'f + {y T y'Y 



In our case of slight roughness, it is sufficient to find only the first three terms of the expansion of the stream 
function ip, Eq. ( A9), in powers of the small parameter e, ip = ipo + eV'i + (2 ^p2 + • • ■ Since all the terms in the operator 
V contain e, the only part of i\) without e is the first term in Eq. (A13) for ipi n h, 



2tt 

V'o (K, y) = (k x ) [exp (iXy) - 1] 



(A15) 



The first order term in ip contains the remaining part of ipi n h and the first order term in the integral ( |AE| ) with 

Vi (k X) d y ) = -£ (k x ) [k 2 x X 2 + k x ] e lXy 

Integration gives 

iX 



ipi (k x ,y) = £ (k x ) 



J^y 



S2 - Si 



{e siy 



(A16) 



The calculation of the second order term requires straightforward integration for much more cumbersome expres- 
sions. However, the general expression for ip2 is not required for further calculations; it is sufficient to have only the 
expression for ip 2 averaged over the surface inhomogeneities, (ip2)- The resulting expression for the stream function 
contains products of the derivatives of the surface profile £' n ) (x) . These products should be averaged over surface 
inhomogeneities using the definition of the correlation function £ (x), Eqs. (Q) 



In the end, after substantial cancellations that accompany the averaging, 



(ip2(k x ,y)) =5(k x ) / dk' x t{K) 



■Sl + 52 

iX 



(iXe iXy + Sl e SlV - s 2 e S2y ) 



(A17) 



(A18) 



Reversing the coordinate transformation of Eq. (A5) and performing the related re-expansion in e, 

'e\ 2 



(v x (K,y-eO)^vW (y)S(k x ) + 



(v^ (k x ,y)) (k x ,y)) + (£) |1««> (v)Hk x ) 



dy 2 



we get for the average velocity 

(v x (x, y)) = exp (iXy) 



00 

/dk 
-^C(k x ){ Sl + S2-iX/2} 



(A19) 



where Si, S2 are given by Eq. (A12). 

The above equations for the stream function allow one to calculate the roughness-driven correction to the dissipation 
of energy and effective friction. 
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Time average of the bulk dissipation per unit area of the wall can be expressed via the stream function ip as 

T)ul 1 



ml 1 

2R A 



dV 



dvj dvk 
dx k dxi 



RA ,dV{Ml« + {* 



where *S? (r,t) = Re [ip (r) e Mt ] , the overline denotes the time average over the period of oscillations, and 

for the statistical average over the random surface inhomogeneities. The time average = ^V'ifeV'ife- 
After the coordinate transformation (|A5|), the attenuation up to the second order term in e reduces to 



(A20) 



stands 



Q = 

Q<°> = 
Q (2) = 
Finally, we get 



ml 

2R 



dy 



Q(°> + e 2 Q (2) 



(0) 



2 -\/2Ay 



2 Re 



A 



The friction force acting on the area unit of the surface is 



F = 2£r, /, 



> n = 



Here n is the unit vector normal to the surface and directed out of the liquid. The square of this force is 



(A21) 



p2 _ f 1 i f 1 2 i TO ~ c S 



J*z f ^ _|_ f 2 <2 

•> Jx ' Jy #2/ 



2 C /2_2 



1 + e 2 f 2 

or, via the stream function, 

In new coordinates (|A5|), this expression reduces to 

</ 2 ) = ((l + 2 e 2 ^)^ (.,0)). 
After separating the real and imaginary parts and expanding in e, we finally get 



(A22) 

(A23) 
(A24) 

(A25) 
(A26) 



</ 2 > 




A 



y/y/kfT^-ki 



(A27) 



Note, that in this problem the friction force introduced by equation ( A23 ) does not determine, after integration 
over the surface, the full energy dissipation. In the case of inhomogeneous rough boundaries there is an additional 
dissipative contribution related to the term with pressure, Prii in the expression for the full force acting on the unit 
area of the surface. If one defines the friction force not via the stress tensor, Eq. (A23), but assumes the experimental 
definition according to F = —Q/uq, then the roughness-driven correction to the friction force will be given by Eq. 
( A22j ) rather than by Eq. ( A27| ). Another anomaly of this problem is that one should always take into account both, 
x and y components of the friction force. 
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FIG. 3: Schematic geometry of the problem. 



APPENDIX B: TWO-LAYER MODEL 



The necessity of using two parameters in the effective boundary condition instead of a single stick-slip length can be 
illustrated by the following simple model. Let us consider tangential oscillations of viscous liquid which is separated 
form a solid substrate by a layer of another liquid with a slightly higher viscosity r\\ > r\ and the same density (see 
Fig. ^|). In effect, we model a rough surface by a layer of viscous liquid with somewhat different properties than in 
the bulk. The model has two parameters: the thickness of the layer d, and dimensionless ratio 7, 

?A Ax 

T = = ^ 1. 

Assuming that the velocity in both liquids is proportional to exp (— iut), we get the following equations of motion: 



d 2 Vl d 2 v 
-iivvi — v\——^ — 0, — vjjv — v—x = 0, 
ay^ ay z 

V! (0) = 1, vi (d) =v(d), 

dvi (d) dv (d) 
Vi J7 = v- 



d y 



dy 



The solution is 



with 



vx (y) = Ae lXl{y - d) + Be- iXx ^- d \ 
v 2 (y) = Ce iX y. 



A,B = Ce 
C 



Xd 



1±7 



-iXd 



cos (Aid) — 17 sin (Aid) 

Time average of the rate of the energy dissipation per unit are consists of contributions from both liquids: 
Q = Q/ + Q11, 



Vi / dy 



Re 



Ce iXd \~ rjxAi 



dvi 
dy 
1 

8V2 



(1 + 7 ) 2 (e^ Ald - l) + (1 - 7) 2 (l - e-^ d ) + 2 ( 7 2 - l) sin (V2A x d) 



Qn = -V / dy 



Re 



\dy e 



= - Ce 



iXd 



2 r/A 
2\/2' 
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If the thickness of the layer d is smaller than the decay length 8, Aid <C 1, 

C - 1 - e «*/*Ad (1 - 7 2 ) , 
|Ce iAd | 2 -> 1 - V2Ad 7 2 + A 2 d 2 7 4 , 



R>e* A " [l-«Ad(l-7 2 



77A 
2Vl 



[1 + A 2 d 2 7 2 (1 - 7 2 



Note, that the condition Aid <C 1 does not necessarily mean that the layer by itself is thin. 
The last two equations show that in this limit 



v (y > d) w Re jitoe 



i(\y— u)t) 



1 



e^/ 4 Ad(l- 7 2 )]}, 



Q 



2\/2 



[l + A 2 dV(l- 7 2 )] 



(Bl) 
(B2) 



Comparison of Eqs.( |5lD , (|_B2| ) with Eqs. (|33|)-(|36|) gives the mapping of the effective viscous layer model onto the 
problem with a rough surface: 

-d(l-7 2 ) = £ eff /R=e 2 £ u 
d 2 7 2 (l- 7 2 ) =e 2 f 2 , 

where £%, £2 are given by the low- frequency equations (p4[) , (psj) . 

In this limit d <C S, the contribution of the layer to the dissipation, Qi, corresponds to the 6— type renormalization 
of the viscosity in the effective boundary condition of Section III with renormalization parameter 







2*i + V2A£ 2 
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